Localization by interference: Square billiard with a magnetic flux 
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Eigenstates and energy levels of a square quantum billiard in a magnetic field, or with an Aharonov- 
Bohm flux line, are found in quasiclassical approximation, that is, for high enough energy. Explicit 
formulas for the energy levels and wavefunctions are found. There are localized states, never before 
noticed in this well studied problem, whose localization is due to phase interference, even though 
there is no or negligible classical effect of the magnetic field. These and related states account almost 
entirely for the magnetic response in certain temperature ranges, and thus have a bearing on the 
experiments of Levy, et. alia 



PACS: 03.65.Sq, 03.65.Ge, 71.70.Di 

The problem of a quantum charged particle subject to 
a magnetic flux and confined in two dimensions, say to a 
square billiard, has attracted much attention. The sim- 
plest fluxes are (1) a uniform magnetic field [UF] and (2) 
an Aharonov-Bohm flux line [ABFL] . One motivation for 
(1) is experimental mesoscopic physics. Levy, et alp have 
made theoretically resistantB measurements of the mag- 
netic susceptibility of an array of such squares. Clearly, 
such a simple problem should be added to the physicist's 
arsenal of known results. 

The Aharonov-Bohm flux lino! has had a rather spec- 
tacular effect on physics. It provides a compelling ex- 
ample of quantum nonlocality. Largely because of this 
example, it is now clearly understood that charged quan- 
tum particles have phase interference effects originating 
from a magnetic field that vanishes in all regions accessi- 
ble to the particle. Although this purely quantum effect 
has no directly corresponding classical physics, it is re- 
lated to the scattering of a classical charge neutral wave 
from a vortexu. Aside from this conceptual motivation, 
an ABFL is a simple idealized extreme case, which like 
the uniform field should be thoroughly understood as 
there are numerous applications. For example, Laugh- 
lin's explanationEl of the quantum Hall effect is based on 
ABFL physics. 

Another reason for interest is that a magnetic flux is 
the most obvious way to break, time reversal symmetry 
of simple systems like billiardsa. The finite flux case is 
in a different universality class from the, zero flux case. 
This has inspired much numerical workQ. Also, square 
symmetry conflicts with the apparent circular symmetry 
of the flux. This, as for the Sinai billiard, suggests chaosS. 
Indeed, the circular billiard in a uniform field is relatively 
trivial and uninteresting. 

Thus, much work has been devoted to these systems. 
It is surprising that there is something new, simple and 
interesting to say about square billiards with a flux. 

We here present some analytic results with numerical 
confirmation. The results are possible because we ex- 
ploit two small parameters. First, we treat the particle 
quasiclassically, that is, its wavelength is short compared 
with the system parameters. Second, we treat the field 
as classically weak, which is automatic for the ABFL. For 
the UF, we assume that e = L/R c is small, where L is 
the square side and R c is the cyclotron radius. 



Nearly all previous work has focused on level statistics, 
or smoothed state densities. In contrast, we find analytic 
wave functions, and simple formulas for their energies. 
The unexpected result is that an important class of these 
eigenstates are localized. Unlike the eigenstates of the 
flux-free case, which uniformly fill the square, some of 
the eigenstates have support in a small fraction of the 
square. The area of this support vanishes at infinite en- 
ergy. Associated with a given extremely localized state 
is a sequence of states which progressively become more 
delocalized, but which share with the localized state a 
simple structure. The most prominent such sequence 
dominates the magnetic susceptibility in the UF case, 
as measured by Levy. 

This localization, like Anderson localization!, is 
caused by phase interference. A bunch of classical phase 
space points mimicking a wave packet in a disordered two 
dimensional region, or in a square, will spread, diffusively 
in the disordered case, linearly in the square, and eventu- 
ally become uniformly distributed. Taking into account 
the phases in impurity scattering, or from the magnetic 
flux, leads to interference effects that suppress the wave 
packet spreading. Many interesting effects of the ABFL 
have been discovered in the last forty years, but we have 
found no previous work that produces a localized quan- 
tum state through destructive ABFL interference. 

In this paper we shall not discuss to its conclusion the 
ideal zero radius ABFL. That limit best establishes quan- 
tum nonlocality, but quantum nonlocality is no longer 
disputed. We have found that the zero radius limit leads 
to interesting but distracting mathematical and numeri- 
cal problems. Basically, there are diffraction effects which 
complicate the results. In this paper we want to avoid 
that issue in favor of understanding the main phase in- 
terference effect. We therefore eventually endow the flux 
line with a finite radius p whose scale we discuss later. 
Initially, however, we suppose p — 0. 

For UF, the cyclotron radius is R c — v/u> c = cp/eB, 
where lo c — eB/mc and p = mv. The momentum p — 
hk = h/X is related to wavenumber k and wavelength A. 
Then e = L/R c — eBLj tick — 2ir<f>i 'fokL. Here (f> is the 
magnetic flux BL 2 , and tfio is the flux quantum hc/e. This 
expression carries over to the ABFL if we identify (f> with 
the flux in the flux line. For the ABFL we may as well 
take \(/>/<f>a\ < 0.5, but for the UF <f> can be large as long 
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as e is small. We choose units such that B = 2n<j)/4>o an( i 
L, Ti and 2m are unity, so that e = B/k. We also suppose 
k » 1, justifyiru^the quasiclassical approximation. 

Our approachEJ utilizes the quasiclassical surface of 
section [SS] method of BogomolnjO. His operator 
T(x,x';E) takes the electron crossing the SS at point 
x' to its next crossing at x, all at energy E — k 2 . The 
SS can be chosen in many ways. To simplify T, we use 
a method of images. Namely, we consider, instead of a 
unit square, x,y € [0,1] <X> [0,1], an infinite channel of 
width 2 obtained by reflecting the original square first 
about x — and then about y — 1 , and finally repeating 
the resulting 2x2 square periodically to |a;| = oo. The 
flux changes sign in neighboring squares. The solutions 
to the square are found from the channel solutions by 
using symmetry, e.g., the solution odd under y 2 — y 
corresponds to Dirichlet conditions at y = 1. 

The SS is taken aSpthe axis y — which is identified 
with y = 2. This givefO 



T(x, x ; E) = 



1 

2nd 



d 2 S(x,x';E) 



dxdx 1 



exp (iS(x, x'; E)) . 

(1) 



Here S — J , p • dr is the action integral along the clas- 
sical path from (x 1 , 0) to (x, 2). Because the field is 
classically weak, this path is approximated by a straight 
line. We immediately find 



S{x, x 1 ) = k\jA+ (x-x'f + $(x, x'), 



(2) 



the flux free result plus $ = (e/c) J A-dr. 

Periodic orbits on the square correspond to straight 
line orbits in the channel from (x 1 , 0) to (x = x' + 2pj 1 q, 
2). Here q is a positive integer and p is a positive or 
negative integer relatively prime to q. Negative and posi- 
tive p are not equivalent if there is a magnetic flux. Such 
orbits correspond to a (p, q) classical resonance which is 
strongly affected by a perturbation. 

Our scheme fuads solutions of Tip = ip by a per- 
turbation theory£2L We first solve J T(x, x')ip(x') = 
e K0 ^ k 'if)(x), treating k as a parameter, and then find the 
energies by solving u>(k) — 2im. Given ip(x), a quadra- 
ture yields the full wave function ^(x,?/). Given 'J, for 
billiards ip (x) oc /dnoc (x,y)/dy\ y= o- 

Specializing to the (1,1) resonance we look for a so- 
lution ip(x) = e lKX u m (x — |), where k — fccos45°, and 
u m varies much more slowly than the exponential. The 
rapidly varying phases in the integral J Tip are station- 
ary at x' — x — 2. This corresponds to diamond shaped 
periodic orbits of the original square whose sides make 
angles of 45° with the x axis. It suffices!!] to evaluate 
$>(x,x') at $>(x,x - 2) = $>(x + 2,x). This is equivalent 
to integrating the vector potential about the closed dia- 
mond loop, and so is independent of gauge. The result 
is that u m satisfies the Schrodinger equation 



where V(x - ±) = -k$(x + 2,x)/C, and C = is the 
length of the diamond orbits. Thus we convert the phase 
$ to a 'potential' V. For UF, 



V(x) — —Bk{\ — 2x 2 )/C; x€[-±,±] 
V(x) 
V(x) 



+Bk(±-2{x + lf)/C- i6[- 
V(x + 2). 
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(4) 



For the (-1,1) resonance, whose orbits are time re- 
versed (1,1) orbits, V(x) changes sign. This would not 
be true if V had its origin in a time reversal invariant 
perturbation of the square, for example, a small change 
of shape. Note also that this periodic extension of the 
x coordinate is similar to use of a 'angle' variable, with 
positive x- velocity v x for x € [0,1], and negative v x for 

ze[l,2]. 

For ABFL, we must choose where to put the flux ling. 
The square center offers some numerical convenienceQ, 
but the results are not very interesting. We place the line 
at (\ : a) which still gives some symmetry and is readily 
compared with the UF case. Then 



V(x) = -Bk/L] x e [-a, a], 

V(x) = +Bk/£; x e [-1 - a, -1 + a] 



(5) 



and for other points in [— |, |], V(x) = 0. This 'square 
well' potential is also continued periodically. Strictly 
speaking, this potential, because of the steps, varies too 
rapidly for the theory to apply. We will modify it below. 

It's clear that Bk is the important parameter in the 
UF case, while for the ABFL, both Bk and a are impor- 
tant. For sufficiently large Bk, Eq. (||), in the UF case, 
will have low lying tight binding harmonic oscillator type 
solutions centered at x = (if B > 0) , with energies 



■hBk/C +(m + ±)y/8Bk/£. 



(6) 



The lowest wave function is approximately uq(x) = 
e -V 'Bk/2Cx wmc ^ j s arbitrarily narrow at large energy. 
For the ABFL 

u m -i(x) w sin (mir(x + a) /(2a)) ; |ac| < a 

w 0; |x| > a (7) 

and the energy is 

E m w -Bk/L + (to + 1) V/4a 2 . 



(8) 



Here the width of the state u is determined by a, and 
the tight binding approximation will be good provided 
Ba 2 k/£ >> 1. The localization of these states in x leads 
to the localization of the full state in the square. 
The total energy is given approximately byEEl 



k n ,m = 2im/ L+E m /k. 



(9) 



+ V(x)u m = E n 



(3) 



Eq. (U) should be solved iteratively. For example, the 
first approximation replaces the k dependence of the 
term E m /k by 2nn/ C. Equivalently, the energy E n ^ n = 
4TT 2 n 2 /£ 2 + 2E m . 



2 



We next give an expression for ^(x,y), shifting the 
origin to the square center: 

*nm(x, V) = (j2 i~ rS K*\ e^ n ^u m (x ~y-\). 

(10) 

Here 1Z : {x,y) — ► (y,—x) is the rotation by 90°. The 
UF Hamiltonian can be taken to be invariant under 1Z, 
while the approximations made for the ABFL induces 
this invariance. Here r labels the rotational symmetry of 
the wave function, i.e. 1Z = i r . Another symmetry gives 
Um(x) = (-l) m M m (-i). The symmetry under transla- 
tion is u(x + 2) = (— l) r u(x). The relation is 

r = — nmod4 + 2(1 — mmod2). (11) 

Thus Eq. (||) holds for all symmetries and successive val- 
ues of n cycle through the representations of 7Z. 

It is easily checked that the & nm vanishes on the 
square sides. This wave function can also be ob- 
tained as a kind of Born-Oppenheimer or channelling 
approximatioi£3. To see that ip m ~ d^! /dn, just observe 
that the derivative, say d/dy at y = — |, need be applied 
only to the rapidly varying exponential. 

Fig. 1. shows numerically calculated | W (x, y) | , which 
is gauge invariant, for n — 62, k rs 2tt62/C rs 138, m = 0, 
B = 25, \jBk w 59. Current density for the same state 
is shown in Fig. 2. For such a well localized u m , each 
term in Eq.(|l0|) dominates one side of the diamond orbit. 
In this case \ip m \ ~ u m . There is interference near the 
square's edges, so the first maximum of 1^1 (e.g. near 
(0,±i)) is about twice as large as |$(±^,±^)|. The 
current density is thus largest close to the middle of the 
square's edges. These states are paramagnetic, that is, 
their current circulates in the opposite sense from that 
of a free particle in the field. 

There is no localized state in the time reversed sense. 
However, the states with m of order ^Bk/C 3 , corre- 
sponding to E m ~ +Bk/2£, energetically at the top of 
'potential' V(x), are diamagnetic. This actually corre- 
sponds to the (1,-1) resonance. These states are not spa- 
tially localized, although they do have a sort of localiza- 
tion in momentum space, as we show in the next para- 
graph. More generally, as m increases, the states become 
more delocalized, and eventually become independent of 
B. This means that for larger m all four terms in Eq. 
( |To| ) make comparable contributions at an arbitrary typ- 
ical point x,y, whereas in the localized case, only one or 
two terms contribute. This gives interference oscillations 
in \ip m \ near \x\ ~ ^ as shown in Fig. 3b). 

We have assumed that n >> m, and that u m is 
slowly varying compared with e lvnx / 2 ^ Expanding, u m = 

Umje 171 ^, where 21 is an integer satisfying (— I) 21 = 
(— l) r . Also u m .i = (— l) m u m ,-i. The unperturbed states 
can be labelled p, q with unperturbed energies (p 2 + q 2 ), 
dropping the factor ir 2 . For r even, Eq. ( |lO| ) is a super- 
position of unperturbed states with quantum numbers 
p = in + I, q = \n — I, where I << ^n. For r odd, 
p = n' + 1' + 1, q = n' — V \ n' ,1' integer. The even r states 



have unperturbed energies \v? + 2l 2 and so are nearly 
'degenerate' to the base energy e„ = \n 2 . In particular, 
they are closer to e n than to the base energy of the next 
representation, e n ±i ~ e„±n. If the perturbation is sym- 
metric under rotation, the next base energies coupled are 
e„^ 4 « e„ ± An. There are, however, many unperturbed 
states with p 2 + q 2 e n . For example, 7 2 + 49 2 = e 70 . 
However, the matrix elements Ti. pq y q i of a smoothly per- 
turbed Hamiltonian, in the unperturbed basis, are small 
if \p — p'\ or \q — q'\ is large. Thus, an interpretation of 
our method, which yields the states of Eq. jlO|), is that 
we effectively diagonalize the Hamiltonian in a basis re- 
stricted to the unperturbed states 'degenerate' with e n 
and close to ^n, \n. This is the case for the UF, and 
indeed, we achieve agreement between full numerical di- 
agonalization, diagonalization restricted to 'degenerate' 
states, and the procedure using the solution of the differ- 
ential equation Eq. (^). 




FIG. 1. Absolute value of a localized wave function on the 
unit square. Uniform field, B=25, n=62, m=0. 

For the ideal, zero radius ABFL, there are significant 
deviations from this scenario. Indeed, most matrix ele- 
ments of the ABFL perturbed Hamiltonian in the unper- 
turbed basis are infinite. However, it is a weak, logarith- 
mic infinity, and our theory seems to capture the main 
shape of the wave function, although at the relatively 
low energies, for which numerical results are available, 
there are significant corrections. We consider these to 
be diffraction corrections, arising from a characteristic 
length shorter than the wavelength. 

We therefore use a finite radius, p, for the flux line. 
The field inside the flux tube is Bo — (j>/Trp 2 . The typical 
angular deflection suffered by a particle traversing this 
field is 59 ~ (<p/ <f>o) / kp. To avoid diffraction we require 
89 to be small. In the numerical work shown, we take 
(f)/(j) = 0.1, and p = 0.01, while k > 140. This hardly 
changes the effective potential of Eq. (||) . An alternative 
and equivalent condition-is to insist that, on the appro- 
priately defined averagec 2 ], the terms in the Hamiltonian 
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satisfy ((eA/c) 2 /2m) << (ep ■ A/mc). 

Fig. 3 is for the ABFL case, with n = 58, m = 0, r = 
and a — j. Fig. 3a) shows uo(x) and uq(— x — 1), its ex- 
tension into x < — J;, reflected. For these parameters, u$ 
is not extremely localized, and extends significantly out- 
side [— i, i]. The remaining plots give \ip m \ — \d^f/dn\ . 
Fig. 3b) plots Eq. ( [l0|), 3c) is from diagonalization in 
the limited basis of the 'degenerate' states; and 3d) is 
from numerical diagonalization in the complete unper- 
turbed basis. To give the flavor of the diffraction effects, 
we show two numerical results for \i/j m \ for a zera_radius 
flux line obtained by a special numerical methodlj: Fig. 
3e) uses the above parameters, and 3f) is for n = 70. 

We finally argue that the (1,1) states dominate the 
orbital susceptibility in a parameter range appropriate to 
experiments. The orbital susceptibility x of a system of 
noninteracting spinless electrons is given by x = dM. /OB 
where the magnetization M. = —d&l(T, /x, B)/dB is that 
of the grand canonical ensemble, 

M(T, f ,,B) = ^2-^f D (E n<m (B)). (12) 

n,m 



Here, uj r — n(2r + l)fcsT and a m — dE m /dB. As an 
example, take fc^T ten times the level spacing J -1 of all 
levels, i.e. k B T = 40tt. Then, u> C/2 « 560. If k F « 500, 
so that the square contains about 4000 electrons, the 
exponential suppression will not be too serious Jar r = 0, 
s = l. Doing the m sum is possible but harderEj. 

Eq. (|lj) is obtained from explicit energy levels rather 
than the action of periodic orbits as found from the 
trace formula, but the result is the same. Integrable sys- 
tems have regularities in their energy levels which lead 
to larger effects in quantities like M as compared with 
chaotic systems. For the square, the (1,1) states have 
the smallest C and also the largest a m . The (2,1) res- 
onance does not couple to a constant field. The (3,1) 
resonances have length £ 3 i = 2\/l0 = VbCn and a 
potential V^i(x) = Vn(x)/3\/5. An interesting case is 
the (1,0) resonance which has Ciq = 2, but this reso- 
nance corresponds to classical orbits enclosing no flux. 
At stronger fields, the curvature of the orbits must be 
allowed for, and eventually this resonance dominateso. 
At very strong fields, such that R c < L, the standard de 
Haas- van Alphen susceptibility is recovered. 



where fo is the Fermi-Dirac distribution function. [The 
canonical ensemble is required in averaging over many 
such billiards, but for a single billiard the grand canon- 
ical suffices!] The chemical potential fi = k 2 F is nearly 
independent of B, since the many states not depending 
much on field act as a heat bath. 
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FIG. 2. Current density of the wave function shown in Fig. 
1. The stream lines are plotted, their density proportional to 
the current density. The sense of the current is paramagnetic. 

It is easiest, no doubt, to estimate Eja. ( ^2[ ) using 
the perturbed Berry- Tabor trace formulaB. However, it 
can be shown that using Eq. (||) in Eq. ( |l2] ) above is 
equivalent to keeping the (1, ±1) periodic orbits and their 
repetitions in the trace formula. Using the Poisson sum 
formula, replace the sum on n in Eq. ( |l2| ) by an integral 
over fc, and do the integral to obtain 

ka 2a m Cuj ( uj r s£\ . £ m „ 

M= £ -^ exp (-wJ sln(£s(fcir ~ 

r,m,s=0 f \ f / r 

(13) 
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FIG. 3. ABFL case. a)-f) u , and b)-f) \d^/dn\, as a 
function of x, for several cases described in the text. 

In conclusion, a new techniqueQ allows analytic solu- 
tion of problems of integrable systems subjected to mag- 
netic flux. Solvable problems include billiards which are 
nearly rectangular, hexagonal or elliptical, as well as cou- 
pled nonlinear oscillators, with uniform or nonuniform 
magnetic flux. Some of the solutions are striking, unex- 
pected and experimentally relevant. Supported in part 
by NSF DMR 9624559 and U. S.-Israel BSF 95-00067-2. 
REP thanks ITP Santa Barbara for support and hospi- 
tality during the early phases of the work. 
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